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Abstract 

We study the effect of a (3+l)-dimensional Chern-Simons electrodynamics on the equations 
governing the dynamics of magnetized plasma and fields. In this model, the Chern-Simons (CS) 
part consists of a dynamical pseudo-scalar field whose space-time derivatives couple with the elec- 
tromagnetic field. We explore the CS corrections to the evolution equation for the magnetic field 
in a plasma with non-zero electrical resistivity. We revisit Cowling's theorem in this context and 
observe that the CS corrections lead to possibly small but non-zero source terms for axisymmetric 
magnetic field. The scalar product of electric and magnetic fields play the role of source of the 
pseudo-scalar field, and therefore, pulsars are likely astrophysical candidates to generate propa- 
gating pseudo-scalar waves. Although aligned electric field gets shorted out by flowing charges in 
large parts of the magnetosphere, there are vacuum gaps in the vicinity of pulsars where strong 
E.B is expected to be present. We derive a wave solution for the pseudo-scalar field generated by 
the time-varying E.B associated with a pulsar. 
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I. INTRODUCTION 

Electromagnetic fields have a dominating influence on plasma behaviour, and therefore, 
any modification of Maxwell's equations is likely to have an interesting observable effect on 
magnetohydrodynamics. Chern-Simons theory, in which the Lagrangian density contains 
the 4-vector potential explicitly, with its non-trivial topological aspects has found many 
applications in physics. It is well known that the (2+l)-dimensional Cherns-Simons (CS) 
theory plays a crucial role in planar physics, particularly in the domain of fractional quantum 
Hall effect and anyonic superconductivity [1,2]. Further more, Witten (1989) established a 
connection between the mathematical theory of knot and link invariants and CS gauge theory 
[3]. With the paper of Carroll, Field and Jackiw (1990) (CFJ90) it became clear that not 
only one can formulate a (3+l)-dimensional CS electromagnetic theory in terms of a fixed 
4-vector p M one can also put constraint on it using radio-data of distant radio-galaxies [4]. 

Thereafter, Jackiw and Pi (2003) (JP03) extended this analysis to construct a (3+1)- 
dimensional CS- modified theory of gravitation, involving again a fixed vector [5]. Setting 
Vii &(j Qx* ^ m a theory based on JP03, it was shown that CS term amplifies the gravito- 
magnetic part of gravity in such a way that speed of a test particle moving in a circular 
orbit around a blackhole reaches saturation at large distances, implying the possibility of 
explaining flat rotation curves observed in galaxies without invoking dark matter [6]. The 
CS term in gravitation also induces a shift in the quantum phase in neutron interferometry 
although the CS correction to phase shift is too small to be observed in current experiments 
[7]- 

To account for the observed late-time acceleration of the universe, it was shown that a 
dynamical scalar-density of weight 1 can be a source of dark-energy that is responsible for 
the acceleration, and that this model entails a natural CS extension of electrodynamics as 
well as gravitation in the (3+l)-dimension, along the lines of CFJ90 and JP03 [8]. In flat 
space-time, the scalar-density appears as a pseudo-scalar field ip that changes sign under 
space-inversion and under time-reversal. 

In the present study, we explore the possibility of Chern-Simons modification of MHD 
equations, and then show explicitly that CS electrodynamics allows one to estimate the 
oscillation characterstics of the pseudo-scalar wave, generated by time-dependent E.B in 
the vacuum gaps of a pulsar magnetosphere, in which the pulsar is taken to be an oblique 
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magnetic dipole rotator. In the following section, we begin with the description of a model 
that incorporates CS modification of Maxwell's electromagnetism. 

II. CHERN-SIMONS MODIFIED MAGNETOHYDRO DYNAMICS 

Two important tensors in physics remain invariant under proper Lorentz transformations 
- the Minkowski metric r] pv and the totally antisymmetric Levi-Civita pseudo-tensor e^ vpc . 
The latter, however, changes sign under space inversion as well as under time reversal. The 
object w pvpa = ip(x)e pvpa is a completely antisymmetric fourth rank tensor under both 
proper and improper Lorentz transformations provided ip is a pseudo-scalar. By supposing 
that w^ypa or, equivalently, tp to be a dynamical field, we had investigated its role as a dark 
energy candidate to explain the observed acceleration of the expansion of the universe [8], 
and the action that was used for ip reduces in flat space-time to, 

where A and B are parameters having the physical dimensions M/T and M/(TL 2 ), respec- 
tively, while ip is dimensionless. 

The pseudo-scalar tp can also be used in conjunction with an approach described in CFJ90 
to write down a gauge invariant Chern-Simons action, 

Scs = J/ ^^A a ^ d*x (2) 

where J has the dimension of T/L. If one assumes pp = 4JcJ^ to be a constant 4- vector, the 
action given by eq.(2) reduces to that described in CFJ90. From the observed correlation 
between polarization and position angles of radio-galaxies, CFJ90 obtained a constraint 
m < 1.2 x 1CT 42 GeV (if one takes the Hubble constant to be H w 72 km/sec/Mpc), where 
m is the mass associated with the Compton wavelength (p^p^) -1 / 2 . As ip is dynamical in 
the present context, this constraint has no bearing on our analysis. However, we expect that 
the magnitude of ip does not deviate too much from unity, and that the dimensionless CS 
coupling constant Jc is small. 

The full action for the CS modified electromagnetism is given by, 

S = / F^F pu d^x-^J fA^x + S^ + S cs , (3) 
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where j M is the electric 4-current density associated with standard matter. 

By extremizing S with respect to and ip, respectively, we obtain the following dynam- 
ical equations, 



J 

c 



8ttc J e^F, 



1 dip 

ip dxl 3 



and, 



1 d 2 

-— - V 2 

c 2 dt 2 



2i/j ' dx» dx u 2 A AA 
In terms of electric and magnetic fields, the above equations take the following forms, 



V.E = 4np e + 16ttc Vg.B 



(4a) 
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where, 



i;(f,t) = (f(r,t)) 2 , 
<7(f,t) = JlogV = 2Jlog/, 



(46) 
(4c) 

(4d) 

(4c) 

(5a) 
(56) 



and p e is the electric charge density. 

For any magnetized plasma, where the fluid bulk speed u is non-relativistic, the magnitude 
of -J^ is only of the order ^ times spatial derivatives of magnetic field, so that it can be 
neglected in comparison to the other terms in order to arrive at, 



„ 4ira ,-> u 

V x B = (E + - x B) - 16vrc 

c c 



1 dg - „ - 
-JLB + VgxE 

C OT 



(6) 



after making use of Ohm's law, with a being the conductivity. Taking the curl of eq.(6), we 
obtain, 



dB 
~dt 



+ V x (B x u) = -V x (77V xB)-Vx 



16™! -^-B + Vg x £ 

C OT 



(7a) 



where 77 = ^ is the plasma electrical resistivity. 
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As g and its derivatives are expected to be small, while the conductivity of cosmic plasma 
is in general large, we may neglect 0(g 2 ), O(gu/o) and 0(g/a) terms, with the result that, 



dB -> -> 

— + V x (B x u) = 7]V 2 B - 16tt?7 

L/6 



4:%dg^ ^fldg\ -; „ _ ' 



(76) 



While going from eq.(7a) to eq.(7b) we have assumed the conductivity a to be a constant 
and have used eqs.(4a)-(4c). 

We note that even after the above approximations, eq.(7b) still implies d ^ = upto 
0(g), for any neutral plasma, since p e = and thereby V.j = 0. Therefore, dynamical evo- 
lution of magnetic field according to eq.(7b) does not result in spurious creation of magnetic 
monopoles, when p e = 0. 

It is interesting to revisit Cowling's theorem [9] in the context of Chern-Simons modified 
electrodynamics. Extending Shu's approach to Cowling's theorem [10] to study the case 
of an axisymmetric distribution of plasma matter and fields, we express the magnetic field 
B(s,z) in circular cylindrical coordinate system (s,(p,z) as, 

— * 

B = B s e s + B^e^ + B z e z , 

where all the components of B are 0- independent. 
With p e = 0, eq.(7b) can be expressed as, 

DB + Biy.u) - (B.V)u = r]V 2 B - lQnr]V x (-^-b) , (8) 



Dt v ' v ' ' ' \cdt 

where, 

D 9 

is the comoving derivative along the flow. 
Using, 



Dt\s 



and, 



which follows from mass conservation, p being the mass density of the plasma, the evolution 
equation (eq.(8)) for the component of the magnetic field takes the following form, 
D (Bs\ 1 r f d (\ d , „ A d 2 B (t> ] („ fldg ~\\ „ du* „ fdu, 



s 
(9) 



First term on the right hand side of eq.(9) leads to the usual decrease in due to 
resistive diffusion. The last two terms consisting of poloidal components of magnetic field 
would also eventually decay [10], and therefore, cannot be an eternal source for B^. The 
second term on the RHS of eq.(9) due to CS correction does raise the possibility of a source 
or a sink (depending on its sign) for the toroidal component, but in normal astrophysical 
situations, its magnitude is expected to be too small to affect the conclusion of Cowling's 
theorem on the impossibility of a steady axisymmetric magnetic field, supported by current 
distributions. However, eq.(7b) could possibly play a significant role in magnetogenesis in 

— * 

the early universe in tiny pockets where fluctuating p e and Vg cause birth of seed B field 
(also see [11]). 

But first of all, is it at all possible to generate time- dependent g? We take this up in the 
next section. 

III. PULSAR AS A SOURCE OF PROPAGATING PSEUDO-SCALAR WAVE 

The evidence for strongest magnetic fields in the universe comes from observations of 
pulsars and magnetars [12, 13]. Pulsars are rapidly spinning magnetized neutron stars with 
high values of angular speed Q and with polar magnetic fields often exceeding 10 12 G. The 
associated electric fields near the surface reach large values causing ions and electrons to be 
pulled out of the surface of the neutron star [14, 15]. A substantial portion of the stripped 
ionized matter that are anchored to the closed magnetic field lines, within the light-cylinder 
of radius c/fi, form a magnetosphere that largely corotates with the pulsar. Rest of the 
plasma matter escape along open field lines as pulsar wind. Since E.B acts as the source 
of the pseudo-scalar field (eq.(4e)), one needs to investigate the possibility of generating ijj 
waves from pulsars as well as magnetars (with \B\ > 10 14 G). 

In this paper, we adopt the oblique rotator model in which the rotating neutron star is 
associated with a magnetic dipole moment jl(t) that maintains an angle a with respect to 
the axis of rotation, as it rotates about it. We assume a geometry in which the spin axis of 
the pulsar is along the z-axis so that its angular velocity is flk. Then, the magnetic moment 
is given by, 

jl(t) = -£> p a 3 [sin a(cos(Qt)i + sm(flt)j) + cosak], (10) 
where a and B p are the neutron star radius and the magnitude of the magnetic fields at the 
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poles, respectively. 

For the rotating magnetic dipole moment, the electric and magnetic fields as measured 
by a stationary observer at the position vector f outside the neutron star but within the 
light-cylinder are given by [16,17], 

E = — %(k x fl(t)) x n + ^-[(k.fl(t))n + (k.n)fl(t) + (fl(t).n)k - 5(fl(t).n)(Ln)n\ 
cr 2 cr 4 

B = ±[3(il(t).h)h-mi (ii) 

where the centre of the pulsar is taken to be the origin of the coordinate system and n is 
the unit vector along r. In spherical polar coordinates, 

fl(t).n = ^B p a 3 K(a,e,<f ) ) (12) 

where, 

K(a, 6, 0) = [sin a sin 6 cos(f2t — 0) + cos a cos 9} (13) 
Then, from eqs.(ll)-(13), 



E.B = 



VLBla? 



p 



a 2 \ 4a 2 



(K cos a — cos 6) ( 1 J H —K cos 6 



(14) 



4cr 5 

In most parts of the magnetosphere, the electric field along the direction of magnetic field 
gets shorted out by the conducting plasma. However, the vacuum gap where electrons and 
positrons are accelerated to very high energies by the aligned electric fields are the sites 
with non-zero E.B [18]. Therefore, the rotating polar and the outer gap regions of the 
magnetoshere, by virtue of eqs.(4e) and (14), are likely to play a significant role for the 
generation of propagating ip field. 

It was shown in [8] that if ip acts as the source of dark energy (DE) that leads to the 
observed late-time acceleration of the universe, then it is necessary that the value of the 



parameter y— ^ is about ^V^D-E ~ 10 28 cm \ ^de being the present dark energy 
density parameter. Thus, for stellar scales which are less compared to the Hubble radius 



« 10 28 cm by a factor of more than fifteen orders, we can neglect the term yf—^f in eq.(4e) 
so that it simplifies to, 



,2 



1 8' 
-— - V 2 

c 2 dt 2 



As / is expected to be close to unity, it can be set to, 

f(r,t) = l + h(r,t) (16) 



where the magnitude of h is much less than one. Then, upto O(Jh), eq.(15) assumes a 
simpler form, 
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1 8' 
-— - V 2 

c 2 dt 2 



h = —E.B, (17) 



which can be integrated to yield the retarded solution, 

v™--^/ ™^:*™ " (18) 

The contribution to the volume integral comes only from the vacuum gap regions where 

— * — * 

E.B is non-zero. There is considerable uncertainty in the literature about the structure of 
vacuum gaps, and the subject of 3-dimensional geometry of polar and outer gap regions is 
an active area of research (see, for example, [19,20]). However, we can infer the nature of 
space-time dependent oscillations of h at distances larger than the size of the magnetosphere 
by substituting, 

(ft \ 
1 - —J (19) 

in eq.(18) and then integrating over r", assuming that E.B in vacuum gaps is given by 
eq.(14). The resulting expression for h has the following form, 
T VtB 2 a & 

h(f, t) « -— — Z—lK^a) + K 2 (a) cos(Q(t - r/c) + 00 + K 3 (a) cos(2fi(t - r/c) + 2 )] 



16nA cr 



(20) 



where, 



Ki(a) = ai + a 2 sin 2 a, 

K 2 (a) = 03 sin a cos a, 

K 3 (a) = a 4 sin 2 a, (21) 

where the constants <f>i, 02, ai, 0,2, ^3 and encapsulate the integrals over the unknown 
gap geometry, and hence not only depend on the exact 3-dimensional structure of vacuum 
gaps but also on n. The first term in the RHS of eq.(20) represents a static h field falling as 
£ while the remaining two terms constitute dipole and quadrupole radiation, respectively. 
Hence, the time variation of h occurs with frequencies fl and 2fl. An interesting but difficult 
problem is to study the propagation of electromagnetic wave through a region where h is 
given by eq.(20). 

Up to first orders in h, eqs.(5b) and (16) lead us to, 



Hence, using eq.(20) in eq.(22) we may estimate the magnitude of at a distance of hundred 

times the light- cylinder radius of a Crab-like pulsar, 

do J 2 / B \ 2 ( a \ 6 / P \ ~ 2 / r \ 

I — I & 5 x 10 42 cm _1 — - x Geometrical factor x — -7— — - — — — — 

l dt l A \10 12 GJ VIOkm/ V33 msj V10 5 km J 

(23), 

with P = being the spin-period of the pulsar. The unknown geometrical factor in eq.(23) 
arises from the integrations over the vacuum gap region (see eq(21)). Although the numerical 
figure in eq.(23) is impressive, until we have estimates of the geometrical factor as well as 
the constants J and A, it is difficult to ascertain the importance of Chern-Simons correction 
to the evolution equation of the magnetic field (eq.(9)). 



IV. DISCUSSIONS 

In our model, the Chern-Simons (CS) action originates due to a coupling between the 
dynamical fourth rank tensor tp(x)e flva ^, the electromagnetic field and the 4- vector potential. 
We derived the MHD equations with CS corrections and showed that the latter gives rise 
to a source term in the evolution equation for the toroidal component of the magnetic field, 
assuming an axisymmetric distribution of both a flowing plasma as well as magnetic field. 
The CS corrections could possibly lead to spawning of seed magnetic fields in small regions 
in the early universe. From the dynamical equations, we found that scalar product of electric 
and magnetic fields acts as a source of the pseudo-scalar field ip. Pulsars and magnetars 
are likely to be strong sources as they are associated with large values of magnetic and 
electric fields. Although aligned electric field is completely screened in most parts of the 
magnetosphere of a pulsar, there are charge-deficit gaps where E.B is expected to be large. 

— * — # 

We obtained explicit wave solutions for tp generated by time-dependent E.B in the vacuum 
gaps of a pulsar magnetosphere. 

However, it is to be noted that the vacuum gaps are fairly close to the neutron star where 
the curvature of space-time geometry cannot be neglected. Strictly speaking, one should 
solve the curved space-time equivalent of eq.(4e) (see, for example, [8]) after having deter- 
mined the background metric in a self-consistent manner using CS gravity around a rotating 
neutron star. Although, in recent studies, weak field corrections to the metric tensor aris- 
ing because of the CS gravity term have been computed assuming the pseudo-scalar field to 
have only time dependence [21,22], it appears that CS gravity does not admit stationary and 
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axisymmetric solution apart from Schwarzschild, flat space-time and two other unphysical 
solutions no matter how the pseudo-scalar field is chosen [23] . In the absence of a Kerr-like 
solution in CS gravity, it is difficult to provide a self-consistent calculation of propagating ip 
field generated from strong gravity regions of a pulsar. 
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